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ABSTRACT 

We present the first non-local iz >0.2) measurement of the cluster-cluster spatial correlation length, using 
data from the Las Campanas Distant Cluster Survey (LCDCS). We measure the angular correlation function for 
velocity-dispersion limited subsamples of the catalog at estimated redshifts of 0.35 < Zest < 0.575, and derive 
spatial correlation lengths for these clusters via the cosmological Limber equation. The correlation lengths that 
we measure for clusters in the LCDCS are consistent both with local results for the APM cluster catalog and with 
theoretical expectations based upon the Virgo Consortium Hubble Volume simulations and the analytic predic- 
tions. Despite samples containing over 100 clusters, our ability to discriminate between cosmological models is 
limited because of statistical uncertainty. 



1. INTRODUCTION 

The spatial correlation function of galaxy clusters provides 
an important cosmological test, as both the amplitude of the 
correlation function and its dependence upon mean interclus- 
ter separation are determined by the underlying cosmological 
model. In hierarchical models of structure formation, the spa- 
tial correlation length, ro, is predicted to be an increasing func- 
tion of cluster mass, with the exact dependence determined by 
ag (or equivalently fio, using the constraint on o-g - f^o from the 
local cluster mass function) and the power spectrum shape pa- 
rameter, r. Low density and low F models generally predict 
stronger clustering for a given mass and a greater dependence 
of the correlation length upon cluster mass. 

The three-space correlation function of clusters was first 
measured for subsamples of the Abell catalog by Bahcall & 
Soneira (1983) and Klypin & Kopylov (1983). Both groups 
found that the correlation function is well-described by a power 
law, ^(r) = (r/r())~^, and obtained a correlation length ro ~ 
25h-^ Mpc with 7~2. Bahcall & Soneira (1983) also observed 
a strong dependence of correlation strength upon cluster rich- 
ness, which was later quantified by Bahcall (1988) and Bahcall 
& West (1992) as a roughly linear dependence of ro upon dc, 
the mean intercluster separation. Postman, Huchra, & Geller 
(1992) and Peacock & West (1992) confirmed the form of the 
correlation function for the Abell catalog in their larger spec- 
troscopic samples, with both studies obtaining ro~20/i~' Mpc 
for clusters with richness class R> 1 . 

While these correlation lengths have strong implications for 
cosmological models, a key problem with interpretation of the 
Abell results is concern that the observed correlation lengths 
are positively skewed by projection effects and sample inho- 
mogeneities (see Sutherland 1988; Dekel et al. 1989; Efstathiou 
et al. 1992; Peacock & West 1992). Several analyses find that 
^((7,7r) is strongly anisotropic, evidence that these effects are 
significant (Sutherland 1988; Efstathiou et al. 1992; Peacock 
& West 1992). Still, the net impact of these factors is unclear. 



Contrary to the concerns raised by these studies. Miller et al. 
(1999) use an expanded sample of Abell clusters to derive cor- 
relation lengths that are consistent with earlier analyses and ro- 
bust to projection effects, and van Haarlem, Frenk, & White 
(1997) argue that projections are insufficient to account for the 
stronger correlation observed in the Abell catalog as compared 
to the APM catalog (Dalton et al. 1992). 

Fortunately, independent constraints on the correlation func- 
tion have arisen as new catalogs with automated, uniform se- 
lection criteria have become available (e.g., Dalton et al. 1992; 
Nichol et al. 1992; Dalton et al. 1994; Nichol, Briel, & Henry 
1994; Romer et al. 1994; Croft et al. 1997; ColHns et al. 2000; 
Moscardini et al. 2000). Some of the recent optical catalogs, 
such as the APM (Dalton et al. 1992), also probe to lower dc 
than the Abell samples, while the X-ray samples provide greatly 
improved leverage for the most massive (highest dc) clusters. 
While systematic variations persist between samples, the exist- 
ing data are generally consistent with ro slowly increasing with 
dc (however, see Collins et al. 2000). 

To exploit the growth of observational data, there has been 
a corresponding theoretical effort to predict the cluster spatial 
correlation function as a function of mass and epoch. Driven 
by cosmological volume N-body simulations (e.g., Governato 
et al. 1999; Colberg et al. 2000; Moscardini et al. 2000) and the 
development of a well-tested analytic formalism (Mo & White 
1996; Sheth & Tormen 1999; Sheth, Mo, & Tormen 1999), a 
theoretical framework has been established that enables deriva- 
tion of quantitative cosmological constraints from the obser- 
vational data. The mass dependence can be studied using ex- 
isting data sets and is typically best matched by low-density 
models (see Croft et al. 1997; Borgani, Plionis, & Kolokotronis 
1999/?; Collins et al. 2000), although systematic uncertainties 
and observational scatter have precluded precision cosmolog- 
ical constraints. In contrast, the redshift dependence remains 
unconstrained because the data have not existed to test the evo- 
lutionary predictions of these models. 
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In this paper we utilize the Las Campanas Distant Cluster 
Survey (LCDCS) to determine the spatial correlation length at 
zc±0.45. We first measure the angular correlation function for 
a series of subsamples at this epoch and then derive the cor- 
responding ro values via the cosmological Limber inversion 
(Peebles 1980; Efstathiou et al. 1991; Hudon & Lilly 1996). 
The resulting ro values constitute the first measurement at this 
epoch of the dependence of the cluster correlation length upon 
dc, probing mean separations similar to previous local optical 
catalogs. Popular structure formation models predict only a 
small amount of evolution from z=0.45 to the present, as il- 
lustrated in section 5. We test this prediction by comparing our 
results with the local observations. 

2. THE LAS CAMPANAS DISTANT CLUSTER SURVEY 

The recently completed Las Campanas Distant Cluster Sur- 
vey, which contains 1073 candidates, is the largest pubhshed 
catalog of galaxy clusters at z > 0.3 (Gonzalez et al. 2001). 
Clusters are detected in the LCDCS as regions of excess sur- 
face brightness relative to the mean sky level, a technique that 
permits wide-area coverage with a minimal investment of tele- 
scope time. The final statistical catalog covers an effective area 
of 69 square degrees within a 78° x 1.6° strip of the southern 
sky (860 X 24.5 /j"^ Mpc at z=0.5 for r2o=0.3 ACDM). Gon- 
zalez et al. (2001) also provide estimated redshifts, Zest, based 
upon the brightest cluster galaxy (BCG) magnitude-redshift re- 
lation that are accurate to ~15% at Zest = 0.5, and demonstrate 
the existence of a correlation between the peak surface bright- 
ness, S, and velocity dispersion, a. Together these two prop- 
erties enable construction of well-defined subsamples that can 
be compared directly with simulations and observations of the 
local universe. 

3. THE LCDCS ANGULAR CORRELATION FUNCTION 

We wish to measure the two-point angular correlation func- 
tion for both the full LCDCS catalog and for well-defined sub- 
samples at z~0.45 that we shall use to constrain the mass de- 
pendence of the correlation length. Gonzalez et al. (2001) find 
that the velocity dispersion is related to the peak surface bright- 
ness of the cluster detection via log cr cx logS(l +z)^ ', and so 
we employ this relation to define subsamples that are roughly 
velocity dispersion limited (see Figure 1). Error in the red- 
shift dependence of this relation can induce a systematic bias 
in construction of our subsamples; however, this uncertainty is 
sufficiently small as to have negligible impact for the LCDCS 
data within the redshift range of the subsamples. We restrict all 
subsamples to z>0.35 to avoid sample incompleteness at lower 
redshift, while the maximum redshift (z=0.575) is set by sample 
size. Subsamples with larger redshift limits (and hence higher 
surface brightness limits) result in smaller samples (A^<100) 
and correspondingly larger statistical uncertainties. 

To compute the two-point angular correlation function, we 
use the estimator of Landy & Szalay (1993), 

nr{nr-\)DD n,-\DR 

w(0)= 1-1, (1) 

n(n-l) RR n RR 

which Kerscher, Szapudi, & Szalay (2000) have demonstrated 
is more robust than other algorithms. In Equation 1, « is the 
number of clusters, rtr is the number of random points, DD is 
the number of cluster pairs with angular separation 6± 60/2, 
RR is the number of random pairs with the same separation, 
and DR is the number of cluster-random pairs. To construct 
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Fig. 1 . — Extinction-corrected surface brightness distribution of the LCDCS 
clusters used to construct subsamples. The solid lines denote the redshift and 
S(l+z)^ ' thresholds of the subsamples. The dotted curve is the surface bright- 
ness threshold of the full LCDCS catalog for E(B-y)=0.05. Candidates at 
the detection threshold in lower extinction regions will lie slightly below this 
curve. Errors in the estimated redshifts move individual data points along 
tracks similar to this curve. 

the random sample, we first generate a list of 44,000 random, 
unmasked locations within the survey region (2000 locations 
in each of the 22 scan sets that comprise the LCDCS survey). 
Next, we reject locations that fail any of the automated cluster 
selection criteria described in Gonzalez et al. (2001), yielding a 
final random sample of 33,886 positions. 
ModelUng the correlation function as a power law. 



w(0) = A„0'-^ = 



1-7 



(2) 



we use a maximum likelihood approach to determine the best- 
fit values for and 7. Similar to Croft et al. (1997) and Bor- 
gani, Plionis, & Kolokotronis (1999^)), we maximize the likeU- 
hood function for an assumed Poisson probabihty distribution. 



/'"' exp(-/i;) 



(3) 



where n; and Hi= < DD > are the observed and expected num- 
ber of pairs in the interval d^. The final results from the max- 
imum likelihood analysis are not dependent upon the exact 
choice of d6, which can be made almost arbitrarily small. To 
determine the best-fit model parameters we use pairs with sep- 
arations of 2' - 5°, and Monte Carlo simulations are utilized to 
determine the associated 1-cr uncertainties. Specifically, we use 
the best-fit parameters as initial values, generate 1000 random 
reahzations of the cluster-cluster pair distribution, and then use 
the distribution of recovered parameter values to quantify the 
observational uncertainties. The angular correlation function 
for the entire LCDCS catalog and for our lowest redshift sub- 
sample are shown in Figure 2, overlaid with best-fit power law 
models. For the full LCDCS catalog, we obtain 7=-1.78±0.10 
and A^=-l .44 ± 0.04 (610=5 1 ±22"). 
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Table 1 

Angular Correlation Lengths 



z 


<z> 


A'' 




7 


/ 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


All 


0.56 


1073 


-1.44±0.04 


1.78±0.10 


0.29 


-1.18±0.04 


0.35-0.475 


0.42 


178 


-1.26±0.17 


2.15±0.19 


0.14 


-1.13±0.19 


0.35-0.525 


0.46 


158 


-1.35±0.33 


2.30±0.33 


0.16 


-1.20±0.33 


0.35-0.575 


0.50 


115 


-1.42±0.47 


2.51±0.45 


0.19 


-1.24±0.45 


All 


0.56 


1073 


-1.73±0.05 


2.15 


0.29 


-1.47±0.05 


0.35-0.475 


0.42 


178 


-1.26±0.13 


2.15 


0.14 


-1.13±0.13 


0.35-0.525 


0.46 


158 


-1.25±0.14 


2.15 


0.16 


-1.10±0.14 


0.35-0.575 


0.50 


115 


-1.14±0.18 


2.15 


0.19 


-0.96±0.18 



Note — (1) Redshift range used to measure the angular correlation. (2) Number of clus- 
ters. (3) Values of logAc^, are for 6 in degrees. (7) Contamination-corrected amplitude 
of the angular correlation function. The second set of parameters in the table is for fixed 
values of the slope 7. 



Table 1 lists information for the full LCDCS catalog and 
the three subsamples at z~0.45, including the redshift range 
spanned by each subsample, the mean redshift of the subsam- 
ple, the best-fit values of and 7, and the estimated frac- 
tional contamination. We also present best-fit values for A^, 
fixing 7=2.15 — equivalent to the best-fit value for the low- 
est redshift subsample and similar to the best fit value for the 
ROSAT All-Sky Survey 1 Bright Sample (Moscardini et al. 
2000, 7=2.1l!g;i). 

Because the LCDCS candidates are not spectroscopically 
confirmed, we must correct the correlation amplitude A^^ for the 
impact of contamination before this data can be used to derive 
the spatial correlation length tq. If we assume that the contam- 
ination is spatially correlated and can be described by a power 
law with the same slope as the cluster angular correlation func- 
tion (a reasonable approximation because for galaxies, which 
are likely the primary contaminant, 7~ 1.8- 1.9 (e.g. Roche & 
Bales 1999; Cabanac, de Lapparent, & Hickson 2000)), then 
the observed angular correlation function is 

Ujie) = A^e'-' = (Aelus,ei-(1 -/)'+Ac„n,ammatK,n/) e'-\ (4) 

where / is the fractional contamination. For detections in- 
duced by isolated galaxies of the same magnitude as BCG's 
at z ~ 0.35 (and identified as galaxies by the automated iden- 
tification criteria described in Gonzalez et al. 2001), we mea- 
sure that Agai is comparable to A^^, the net clustering ampU- 
tude for all LCDCS candidates at 0.3<z<0.8. For detections 
identified as low surface brightness galaxies (including some 
nearby dwarf galaxies) we measure Alsb — lOA^j,. While these 
systems are strongly clustered, we expect that they comprise 
less than half of the contamination in the LCDCS. For mul- 
tiple sources of contamination the effective clustering amph- 

filde Acontamination = Y.^ifiKYjfif^ effective CluSterfug 

strength of the contamination is Acontamination % 2.5A^^ even in- 
cluding the LSB's. 

The last column in Table 1 gives the contamination-corrected 
value of A^, under the assumption that Acontamination = A^^ . If the 
effective correlation amplitude of the contamination is in the 
range Acontamination = {0-2.5)A^^ (which corrcspouds to uncor- 
related contamination for the lower hmit), then the systematic 
uncertainties in A^^ and the corresponding uncertainties in the 



ro values derived in §5 are < 1%. These uncertainties, which 
are small because the contamination only contributes to the ob- 
served correlation function with weight are far less than the 
statistical uncertainties in all cases. Uncertainty in the frac- 
tional contamination of the catalog yields a larger 6% system- 
atic uncertainty in ro, which is also less than the typical statis- 
tical uncertainty. 

4. THE SPATIAL CORRELATION LENGTH 

The observed angular correlation function can be used to de- 
termine the three-space correlation length if the redshift distri- 
bution of the sample is known. This is accomplished via the 
cosmological Limber inversion (Peebles 1980; Efstathiou et al. 
1991; Hudon & Lilly 1996). For a power-law correlation func- 
tion with redshift dependence /(z), 

^('•)=(^) ^x/(z). (5) 

The corresponding comoving spatial correlation length is 
roiz) = I'ofizy^'', and the Limber equation is 
^ r(7/2) 

J'f(dN/dzfE(z)DA(zy-^f(z)(l +z)dz 

} ^2 ' 

(J^^(dN/dz)dz^ 

where dN/dz is the redshift distribution of objects in the sam- 
ple, D^iz, Oo, Oa) is the angular diameter distance, and 

Eiz) = [0(l+^)3 + 0«(l+^)-'+^^A]l/^ (8) 

as defined in Peebles (1993). Because httle evolution in the 
clustering is expected over the redshift intervals spanned by our 
subsamples (see the Appendix and Figure 3), f{z) can safely be 
pulled out of the integral. 

We estimate the true redshift distribution of clusters in each 
LCDCS subsample, dN/dz, based upon the observed distri- 
bution of estimated redshifts, dNohs/dz- Because the redshift 
distribution of the full LCDCS catalog is slowly varying over 
the redshift interval probed by our subsamples, we are able to 
derive approximate models for the dN/dz of the subsamples by 
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Table 2 
Spatial Correlation Lengths 

z ACDM(i2o=0.3) OCDM (f^o = 0.3) TCDM(r=0.lf 

range dc dc ro dc 

0.35-0.475 38!4 14.7!^-^ 33^8 12.9!H 309 11.8!i| 

0.35-0.525 46.3 Ki^^jl 40.6 14.2!2| 36.9 12.9!|| 

0.35-0.575 58.1 18.2^^-^ 50.8 15.8^^;^ 46.0 14.4^^j 

Note — The units for dc and ro are ft"' Mpc. Listed uncertainties are statistical. 
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Fig. 2. — Angular correlation functions for the entire LCDCS catalog and 
for a volume-limited subsample with 0.35 < z < 0.475. The solid lines are 
best-fit power law models, with the best-fit parameters given in Table 1 . The 
data shown are plotted with bins of width A log 9=0.2; for the actual maximum 
likelihood analysis of the subsamples we use A log 6=0.005. 



convolving the dNohs/dz with Gaussian scatter {a^/z~0.\A at 
z=0.5; Gonzalez et al. 2001), which eliminates the sharp red- 
shift hmits imposed on Zest in construction of the subsamples. 
The rigorously correct method for obtaining the dN/dz is to 
c/econvolve the scatter; however, our approximation is suffi- 
ciently accurate for use in the Limber inversion because of both 
the slow variation of the LCDCS redshift distribution at this 
epoch and the weak dependence of the Limber inversion upon 
dN/dz. To test the vaUdity of this approach, we also try mod- 
eling the dN/dz using the theoretical mass function of Sheth & 
Tormen (1999) convolved with redshift uncertainty. Comparing 
these two methods we find that the derived spatial correlation 
lengths agree to better than 3% for all subsamples. This result 
is insensitive to the exact mass threshold assumed for the theo- 
retical mass function, again due to the weak dependence of the 
Limber inversion upon dN/dz. 

5. RESULTS AND COMPARISON WITH LOCAL DATA 

Table 2 lists the correlation lengths (ro) and mean sepa- 
rations {dc) that we derive for each of the three subsamples, 
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Fig. 3. — Comparison of the LCDCS data with local samples for a ACDM 
cosmology with f2o=0.3. The error bars on the LCDCS data correspond to 
the 1-cr statistical uncertainty. Overlaid are the results from the Virgo Consor- 
tium Hubble Volume simulations (Colberg et al. 2000) and analytic predictions 
(r=0.2) for z=0 and z=0.45. These fiducial curves do not account for the im- 
pact of uncertainty in cluster masses (see Figure 4). 

treating both 7 and A^^ as free parameters. Both the ro and 
dc values are cosmology-dependent, and so we list the de- 
rived values for ACDM (O()=0.3), OCDM (l^o=0.3), and 
rCDM cosmologies. Values of dc in Table 2 are computed 
as dc = n^l^ = [A?(l -/)/y]"'/^ where N is the number of 
clusters, / is the contamination rate, and V is the effective co- 
moving volume for a given subsample. Including the impact 
of Gaussian scatter upon the redshift distribution, the effective 
comoving volume for a subsample with redshift bounds zi and 
Z2 is 

r- dV 

V = dnx / dzp^jp^ (9) 



dfldz 



The effect of the redshift uncertainty is contained in p, which is 
the probability that a cluster at redshift z will be observed to 
have an estimated redshift within the redshift range spanned by 
the subsample. For Gaussian redshift scatter, 

Piz)- 



t p-(.z'-zfl2al 



(10) 



Inclusion of p increases the computed values of dc by -^3%. 
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We compare the LCDCS results with several local studies to 
assess the degree of evolution between z=0.5 and the present 
epoch. We find that the LCDCS ro values are statistically con- 
sistent with results from the Edinburgh-Durham Galaxy Cat- 
alogue (Nichol et al. 1992), APM survey (Croft et al. 1997), 
and the MX Survey (northern sample, Miller et al. 1999), but 
smaller than the correlation lengths found by Peacock & West 
(1992) for their subsample of the Abell catalog. The correla- 
tion lengths from these studies are plotted in Figure 3, where 
the LCDCS data points are shown for an assumed ACDM cos- 
mology (1^0=0.3). The lowest dc data points for the XBAC cata- 
log from Abadi, Lambas, & Muriel (1998), which probe higher 
masses than our study, are also plotted. We also include fiducial 
theoretical curves to illustrate the predicted degree of evolution 
between these epochs. The thick solid curve shows the z=0 re- 
sults from the Virgo Consortium Hubble Volume simulations 
(Qo = 0.3, = 0.7, as = 0.9,r=0.17). The other two curves 
correspond to analytic predictions for the dependence of ro on 
dc (or equivalently, on number density) at the two epochs, using 
the analytic model of Sheth & Tormen (1999) with r=0.2. 

It is evident in Figure 3 that the weak evolution in the corre- 
lation length predicted for this ACDM model is consistent with 
the observations; however, for a more accurate comparison we 
should incorporate into the models the impact of observational 
scatter in the subsample mass thresholds for the different data 
sets. We restrict our attention to the APM and LCDCS cata- 
logs in this analysis. For the APM catalog, we use the corre- 
lation lengths from Croft et al. (1997), taking the velocity dis- 
persion scatter about the sample thresholds to be dloga=0.13 
(based upon the data set of Alonso et al. 1999). We estimate 
the mean redshift of the APM samples to be <z>~0. 1 . For the 
LCDCS catalog we use the surface brightness-velocity disper- 
sion relation and surface brightness scatter from Gonzalez et al. 
(2001). The data and resulting models that include this scatter 
are shown in Figure 4a-c for the three cosmologies discussed 
in this paper. Both the APM and LCDCS data are consistent 
with the low-density models. For the plotted curves, the com- 
bined data sets yield reduced values xj.=0.7 and x^=0.6 for 
the ACDM and OCDM models, respectively.^ In contrast, the 
plotted rCDM model systematically underpredicted the ro val- 
ues for both data sets (x^=4. 1 for r=0.2). rCDM can be made 
to better match the data by decreasing F; however, it is not pos- 
sible to make rCDM simultaneously consistent with the clus- 
ter correlation length data and the galaxy power spectrum (e.g. 
Eisenstein & Zaldarriaga 2000). 

6. SYSTEMATIC UNCERTAINTIES 

It is important to consider whether any systematic biases 
have the potential to qualitatively alter our results. Table 3 
summarizes the systematic effects discussed thus far in the 
text, as well as several additional potentially important factors. 
The most significant potential systematic bias is attributable 
to the impact of large scale structure. The LCDCS serendip- 
itously includes the most X-ray luminous cluster known (RX 
J1347.5-1145 at z=0.45; Schindler et al. 1995). Three of the 
ten most massive LCDCS candidates with estimated redshifts 
Zes;<0.58 have projected separations of less than 50 Mpc 
relative to RX J1347.5-1145, indicating that this region likely 
contains a massive supercluster. To estimate the sensitivity of 
our results to the presence of this region within the survey, we 

^ The bulk of the weight in computing these xi values comes from the APM data 
uncertainties associated with these measurements. 



recompute the correlation lengths excluding candidates within 
lOO/i"' Mpc of RX J1347.5-1 145. For all three subsamples the 
results change by less than 15% (i.e. less than I -a), and for 
the lowest redshift subsample the correlation length actually in- 
creases slightly when this region is excluded. We thus conclude 
it hkely that no individual supercluster should systematically 
alter our results by more that 15%. 

TABLE 3 

Systematic uncertainties 

Issue Impact 

ro 

Large scale structure ^15% 
Fractional contamination 6% 
Clustering strength of contamination ^10% 
Model for dN / dz in Limber equation 3 % 

Extinction-dependent selection effects 1 % 

Redshift dependence of E < 1 % 

Fixing 7. Would increase ro for subsamples 0-10% 
(larger change for higher dc bins) 

dc 

Uncertainty in effective comoving volume <3% 



A second key concern, due to the cluster detection method, is 
whether there exist extinction-dependent selection effects that 
impact the derived correlation lengths. As a test we recompute 
ro for the lowest redshift subsample using only regions with 
E(B-y)<0.06, which reduces the sample from 178 to 89 clus- 
ters. The recovered correlation length for these low-extinction 
regions is 13.2 (vs. 14.7 for all 178 clusters). This 10% change 
is less than the statistical uncertainty, indicating that our results 
are not strongly sensitive to variations in the galactic extinction 
within the survey region. 

Of the factors previously discussed in the text, the greatest 
potential impact upon our results can be achieved if the a priori 
assumption is made that all subsamples have the same power 
law slope, 7. By making this assumption we test whether our 
results are robust to any factor that may be artificially enhanc- 
ing the best-fit values of 7. If we fix 7=2.15, then the derived 
ro values for three subsamples increase by 0%, 6%, and 15% 
(in order of increasing maximum redshift for the subsamples). 
While these changes are less than the statistical uncertainties, 
they do act to steepen the slope of the observed r^-dc relation 
and thus slightly degrade the agreement with low-density mod- 
els. 

We next consider two additional issues that might qualita- 
tively change the results of this paper. First, what if the red- 
shift uncertainty is significantly larger than indicated by Gon- 
zalez et al. (2001)? In this case the redshift ranges spanned 
by our subsamples are larger — and more similar to one an- 
other — than our current best assessment. If true, then we are 
currently underestimating ro. Quantitatively, if (Tj./z=0.25 in- 
stead of 0. 14, then the actual values of ro would be ~2 Mpc 
larger. On the other hand, neither the mean separation dc for a 
given sample or the steepness of the relation between ro and dc 
changes perceptibly. Second, what if the contamination fraction 
is uniform in all the subsamples rather than being an increas- 

t. The LCDCS provides only slight additional leverage due to the large statistical 
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Fig. 4. — Comparison of the LCDCS and APM data sets with analytic predictions. Different panels correspond to the different cosmologies assumed in deriving 
the LCDCS correlation lengths (see Table 2). The analytic models include the impact of uncertainty in clusters masses for the two data sets, which causes these 
curves to be slightly shallower than the fiducial models in Figure 3. 



ing function of redshift? Qualitatively, this would imply that 
part of the observed weak dependence of ro upon dc is simply 
due to our contamination correction. Quantitatively, if we take 
the mean contamination rate to be 16% for all three subsam- 
ples, then the ro values(ACDM) change from (14.7,16.2,18.2) 
to (15.0,16.2,17.7) /i"' Mpc - a minor effect compared with 
other uncertainties. 

From our analysis of the above systematics, we conclude that 
the two largest systematic uncertainties in this work arise from 
large scale structure and potential underestimation of the red- 
shift uncertainty. Large scale structure can impact the derived 
values of ro at the 1 - cr, while underestimation of the redshift 
uncertainty could yield correlation lengths that are systemati- 
cally too small (by up to ~2 h"^ Mpc if a^/z=025). 

7. DISCUSSION AND CONCLUSIONS 

The Las Campanas Distant Cluster Survey is the largest ex- 
isting catalog of clusters at z>0.3, providing a unique sam- 
ple with which to study the properties of the cluster popula- 
tion. We have used the LCDCS to constrain the cluster-cluster 
angular correlation function, providing the first measurements 
for a sample with a mean redshift z>;0.2. From the observed 
angular correlation function, we derive the spatial correlation 



length, ro, as a function of mean separation, dc- We find that the 
LCDCS correlation lengths are in agreement with results from 
local samples, and observe a dependence of ro upon dc that is 
comparable to the results of Croft et al. (1997) for the APM 
catalog. This clustering strength, its dependence on number 
density, and its minimal redshift evolution are consistent with 
analytic expectations for low density models, and with results 
from the ACDM Hubble Volume simulations. Consequently, 
while statistical uncertainty limits our ability to discriminate 
between cosmological models, our results are in concordance 
with the flat ACDM model favored by recent supernovae and 
cosmic microwave background observations (e.g., Riess et al. 
2001; Pryke et al. 2001; de Bernardis et al. 2001). 
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APPENDIX 

ANALYTIC MODELS FOR THE SPATIAL CORRELATION FUNCTION 

Analytic models for the spatial correlation function are based upon the Press-Schechter formalism (PS; Press & Schechter 1974) 
formalism. Mo & White (1996) were the first to use this formahsm to derive the expected correlation for mass-dependent bias, 
finding that if the observed correlation function is defined as 

^{r,z,M) = b\z,mm{r,z), 

where is the average mass correlation function, 

r°° sinfcr 

Ur,z) = Dl(z) Pik)^d% (A2) 



(Al) 



then the bias, b(z,M), is given by the equation 



b(z,M)=l + 



Sc(z) 



kr 
1 



aHM) ScizY 



(A3) 



Modeling of observed correlation functions requires calculation of the effective bias, beffiz,M), which is the average of b(z,M) 
weighted by the mass function of the clusters used to compute the correlation function (Mo & White 1996; Matarrese et al. 1997; 
Baugh, Cole, & Frenk 1998), 

Cb(MMM)dM 

CnmdM ' ^^^^ 

where m is the mass limit of the cluster sample. 

This prescription can also be easily extended to incorporate the Sheth & Tormen modifications to PS by replacing equation (A3) 
with 



biz,M)=l + 



1 



SciO) 



aS^iz) 



1 



2p 
SciO) 



2p- 



(A5) 



where a=0.707 and p=03 (Moscardini et al. 2000). This is the analytic model used in §5. For comparison, we also overlay results 
from the Virgo Consortium Hubble Volume simulations (Colberg et al. 2000) in Figure 3a. 



